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Abstract 

Using RHESSI data, we have analyzed some 172 hard X-ray peaks during 53 
solar flares which exhibited a double-footpoint structure. Fitting both footpoints 
with power-laws, we find that spectral index differences A7 range mostly between 
to 0.6, and only rarely go beyond. Asymmetries between footpoints were not 
observed to be significantly dependent on their mean heliographic position, their 
relative position with respect to each other, nor their orientation with respect 
to the solar equator. Assuming a symmetric acceleration process, it is also clear 
that differences in footpoint spectral indices and footpoint flux ratios can seldom 
be attributed to a difference in column densities between the two legs of a coronal 
loop. Our results corroborate better the magnetic mirror trap scenario. More- 
over, footpoint asymmetries are more marked during times of peak HXR flux 
than when averaging over the whole HXR burst, suggesting that the magnetic 
configuration evolves during individual HXR bursts. We observed also a linear 
correlation between the peak 50-keV fiux and the peak GOES l-SA channel fiux, 
and that HXR burst duration seem correlated with loop length. 
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1. Introduction 

Solar flares are believed to be powered by magnetic reconnection high in the 
corona, which accelerate particles. Particles, and in particular electrons, travel 
down field lines and emit bremsstrahlung hard X-rays (HXR) as they penetrate 
the denser chromosphere. Hence, the HXR radiation from the footpoints con- 
tain much information about the accelerated electrons (though convoluted with 
transport effects), such as energy content (?; ?; ?). 

The first observations of HXR footpoints were made by ? (?) on SMM. Yohkoh 
HXT later characterized the "standard" flare model: two footpoints and an 
above-thc-looptop source (?). The above-the-looptop source is rarely observed, 
although that could be due to observational constraints. 
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Statistical studies of HXR footpoints (and looptop) sources' spectral indices 
using the Yohkoli satellite (?; ?) have been carried out in the past. Some of their 
results seem to indicate that the spectral indices of two neighbouring footpoints 
could differ by as much as 1 or even 2! Such differences cannot be explained 
easily by transport mechanisms (see e.g. Appendix A). Yohkoh HXT's results 
were compromised by the fact that it had only had 4 energy channels, and had 
to deal with issues like thermal contamination of the lowest channel(s), and the 
sometime poor statistics of the upper ones. 

The Ramaty High Energy Solar Spectroscopic Imager (RHESSI, ?, ?) offers an 
unprecedented combination of spectral resolution (1 keV in the 3-100 keV range), 
spatial resolution (2.3"), temporal evolution (w2s), and sensitivity comparable 
to Yohkoh HXT's. Previous RHESSI imaging spectroscopy results include: ? 
(?), ? (?), ? (?). A similar time variation of spectra in footpoints is observed, 
although in some cases, the spectrum in one footpoint is steeper than in the other 
one (by about 0.3 in spectral index, in ? (?). ? (?) suggest that the discrepancy 
could be due to a difference in column densities of the electron population, as 
they propagate down an asymmetric loop (cf Appendix A). 

All RHESSI papers so far discuss only a single event, or a few events, but no 
statistical study of the footpoints has been done so far, that exploit RHESSI's 
large database of observed flares. These results should help constrain energy 
release and particle acceleration in the flare model. 

2. Event Selection 

RHESSI has been in orbit for six years now. Since its launch on 5 February 
2002, it has observed more than 20000 flares. For this study, the strongest flares 
will be taken into considerations, as they have better count statistics. More 
specifically, flares which show substantial emission above 50 keV (where solar 
flare thermal components are always negligible), enough to produce images of 
good enough quality to be used in imaging spectroscopy (2000 counts is about 
the minimum for a reliable single image reconstruction with two sources). To 
achieve our objectives, we will use the simplest possible events, those that show 
only two footpoints. 

The HESSI Experimental Data Center (HEDC, ?, ?) was used to find our 
events. It was queried for all flares between 13 February 2002 and 1 July 2006 
which had 2 or more sources, and with peak GOES flux above Ml.O level. More 
than 1100 flares corresponded to that description. Each were individually ex- 
amined, in particular the time vs. energy panels of RHESSI images that HEDC 
automatically produces for each flare (one minute accumulations over whole flare 
duration). The flares that were retained were those that visually displayed two 
footpoints in HEDC images above 50 keV. 53 flares were kept (see Table 1). 

3. Method of Analysis 

Imaging spectroscopy using CLEAN (?) and the OSPEX spectral analysis soft- 
ware were employed. Imaging was done using collimators 3 to 8, yielding a formal 
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Table 1. Flares studied, GOES X-ray class, and position on the sun 
in both arcseconds from suncenter and heliographic coordinates. 



Flare Date and GOES Location on Sun 

number time class X ["] Y[ ] Lat/Lon 



1 


2002/02/20 


11 


:06; 


:12 


07. 


5 


904.4 


261.5 


N13 W73 


2 


2002/03/17 


10: 


:15; 


:36 


Ml 


.4 


-342.2 


-239.2 


S20 E22 


3 


2002/03/17 


19; 


:28; 


:44 


M4 


.4 


-264.8 


-232.9 


S20 E17 


4 


2002/04/10 


12 


:28; 


:04 


M8 


.8 


-20.8 


434.5 


N20 EOl 


5 


2002/04/10 


19; 


:02; 


:48 


Ml 


.8 


-348.2 


377.5 


N17 E22 


6 


2002/05/31 


00; 


:07; 


:08 


M2 


.4 


-817.5 


-475.4 


S30 E87 


7 


2002/06/01 


03; 


:53; 


:40 


Ml 


.6 


-414.8 


-293.0 


S18 E27 


8 


2002/06/02 


11 


:44; 


:32 


Ml 


.0 


-148.8 


-300.3 


S18 E09 


9 


2002/07/17 


07; 


:02; 


:48 


M9 


.2 


288.5 


246.5 


N19 WIS 


10 


2002/07/18 


03; 


:32; 


:56 


M2 


.5 


421.4 


264.5 


N20 W28 


11 


2002/07/23 


00; 


:28; 


:04 


X5. 


1 


-868.1 


-235.3 


S12 E70 


12 


2002/07/29 


10; 


:39; 


:08 


M5 


.1 


238.9 


-291.3 


S12 W14 


13 


2002/07/31 


01; 


:48; 


:40 


Ml 


.4 


558.4 


-220.2 


SOS W36 


14 


2002/08/03 


19; 


:04; 


:36 


XI. 


2 


899.8 


-265.8 


S15 WSO 


15 


2002/08/21 


01; 


:39; 


:16 


Ml 


.6 


689.2 


-246.7 


SIO W47 


16 


2002/08/22 


01; 


:52; 


:00 


M5 


.9 


798.4 


-266.5 


S12 W59 


17 


2002/09/08 


01; 


:39; 


:08 


Ml 


.6 


-908.8 


-193.1 


S09 E75 


18 


2002/09/27 


03; 


:34; 


:28 


Ml 


.0 


-694.7 


142.9 


N13 E47 


19 


2002/12/04 


22; 


:47; 


:00 


M2 


.5 


-836.7 


227.7 


N13 E61 


20 


2003/04/23 


01; 


:01; 


:56 


M5 


.1 


261.8 


366.8 


N17 W16 


21 


2003/05/29 


01; 


:04; 


:40 


XI. 


2 


494.4 


-106.6 


S07 W31 


22 


2003/06/17 


22 


:53; 


:40 


M6 


.8 


-790.4 


-138.6 


S07 E57 


23 


2003/07/17 


08; 


:19; 


:46 


09. 


7 


-206.9 


178.2 


N15 E13 


24 


2003/10/23 


08; 


:47; 


:20 


X5. 


4 


-904.8 


-317.4 


SIS ESI 


25 


2003/10/24 


02; 


:48; 


:32 


M7 


.7 


-865.1 


-341.1 


S19 E71 


26 


2003/10/29 


20; 


:43; 


:20 


XIO 


90.8 


-381.3 


SIS W05 


27 


2003/11/01 


22 


:33; 


:04 


M3 


.3 


818.6 


-253.8 


S13 W60 


28 


2003/11/03 


09; 


:49; 


:16 


X3. 


9 


917.8 


130.3 


NOS W73 


29 


2003/11/04 


19; 


:33; 


:56 


M9 


.1 


900.8 


-335.0 


S19 WSl 


30 


2004/01/06 


06; 


:22; 


:32 


M2 


.7 


-972.6 


88.3 


NOS ES9 


31 


2004/01/07 


10; 


:22 


:12 


M3 


.7 


-930.2 


117.0 


N05 E73 


32 


2004/04/06 


13; 


:22; 


:48 


M2 


.3 


-261.5 


-170.6 


S16 E16 


33 


2004/07/13 


00; 


:15 


:26 


M6 


.8 


654.7 


181.6 


N14 W45 


34 


2004/07/23 


21 


:19; 


:28 


Ml 


.8 


125.5 


6.6 


N05 W07 


35 


2004/09/12 


00; 


:33; 


:44 


M4 


.8 


-706.2 


-35.9 


N02 E47 


36 


2004/10/30 


03; 


:30; 


:09 


M3 


.5 


316.8 


145.2 


N12 W19 


37 


2004/10/30 


16; 


:24; 


:26 


M6 


.0 


427.6 


139.5 


N12 W26 


38 


2004/10/31 


05; 


:32; 


:03 


M2 


.4 


540.6 


152.6 


N12 W34 


39 


2004/11/03 


03; 


:30; 


:52 


Ml 


.6 


-674.1 


94.7 


NOS E44 


40 


2004/11/06 


00; 


:30; 


:48 


M9 


.5 


-79.3 


83.4 


NOS E04 


41 


2004/11/06 


01; 


:42 


:34 


M3 


.7 


-27.4 


68.1 


N07 EOl 


42 


2004/11/10 


02; 


:09; 


:44 


X2. 


6 


700.1 


91.4 


N07 W46 


43 


2004/12/01 


07; 


:10; 


:16 


Ml 


.2 


-335.8 


128.5 


NOS E20 


44 


2005/01/15 


06; 


:28; 


:31 


M8 


.4 


-106.1 


295.5 


N12 E06 


45 


2005/01/15 


22 


:48; 


:24 


X2. 


7 


103.1 


306.7 


N13 W06 


46 


2005/01/17 


10; 


:00; 


:23 


X3. 


9 


430.0 


292.3 


N13 W26 


47 


2005/01/19 


08; 


:12 


:40 


M8 


.7 


708.9 


283.6 


N13 W4S 


48 


2005/01/19 


10; 


:21 


:08 


M2 


.5 


679.6 


339.9 


N16 W46 


49 


2005/01/20 


06; 


:44; 


:44 


X7. 


1 


818.5 


256.0 


N12 W59 


50 


2005/07/13 


14; 


:14; 


:14 


M2 


.7 


909.4 


168.9 


Nil W78 


51 


2005/08/22 


01; 


:11 


:54 


M2 


.3 


717.8 


-248.7 


SIO WSO 


52 


2005/08/22 


17; 


:07; 


:34 


M5 


.2 


801.4 


-241.6 


Sll WS9 


53 


2005/08/23 


14; 


:46; 


:21 


M2 


.7 


883.8 


-219.0 


Sll W71 



FP.tex; 21/11/2011; 1:27; p. 3 



Saint-Hilaire et al. 



6-12 keV 
50-100 keV (xl5) 



] 000 ='V./^JV\AMAAA/ 




14:20 14:2b 1 

Start Time (13-Jul-05 14:10:00) 



Recorstructed RHESSI Image 



220 I 

^ I 

g 200 I 

o I 

> 180 I 
% 1 60 I 



900 920 940 960 980 
Heliocentric X (orcsec) 



30 40 50 

Energy [keV] 



Figure 1. RHESSI time profile (top), 30-33 keV image (bottom left) and footpoint spectra 
(Bottom right) of the 13 July 2005 M2.7 flare. The gray bars represent the time intervals chosen 
for our analysis: a "Peak Times" 16-second interval around 14:14:10, and two "Whole Peaks" 
from 14:13:47 to 14:19:50, and from 14:21:22 to 14:24:41. The third HXR peak (Kil4:30) was 
discarded as it was faint, and there was an attenuator state change (Lin et al., 2002) during 
it, further complicating the analysis. The sun was eclipsed by the Earth until ^14:13. 



image resolution of 7" FWHM. The time intervals and energy intervals were 
chosen as follow. 

For each flare, two types of time intervals were used: 

• "Peak flux" time intervals: These arc four RHESSI spin periods long (each 
RHESSI spin period being about 4s long), centered at the time of peak 
HXR (above 50 keV) flux. The later is found using RHESSI Observing 
Summary data (?). The peak flux time interval was taken to be this time 
of peak flux plus or minus two RHESSI spin periods (which are sa4s long) . 
Of course, there can be only one such peak flux time interval per flare, 
resulting in 53 such peaks in our study. Taking the flare of 13 July 2005 
as an example (Figure 1), the time interval of accumulation was about 
14:14:05 to 14:14:22. 

• "Whole peak" time intervals: Strong non-thermal peaks appearing in RHESSI| 
spectrograms (or dynamic spectra) were selected over their whole time 
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interval, defined as the time the HXR flux (>50 keV) is greater than 50% 
of its peak value. There can be many such "Whole peak" time intervals 
within the same flare. And of course, for each flare, one of the "Whole 
peak" time interval envelops that flare's "Peak flux" time interval. Again 
taking the flare of 13 July 2005 as an example (Figure 1), the time intervals 
of accumulation were about 14:13:47 to 14:19:50, and 14:21:22 to 14:24:41. 
In a few cases, time intervals did not contain two footpoints (but only one, 
or sometimes three or more), and were hence discarded from the study. 

For each time interval, the energy binning was chosen using the following 
semi-empirical approach: 

• The start (lowest) energy was visually chosen by inspection of the RHESSI 
spectrogram: it is taken to be the point where the non-thermal emission 
starts to be clearly stronger than the thermal component. 

• The energy binning was taken to be pseudo-logarithmic, which each energy 
bin having at least 2000 counts above background, and the bin width being 
between 5% and 20% of the bin value. This was crudely approximated 
using Observing Summary 4 second data rates. The end (highest) energy 
was taken to be when the next energy bin could not achieve 2000 counts 
above background. We had 4 to 18 (typically 10) energy bins to fit and 
obtain spectral indices and fluxes with. 

Finally, for the results presented and discussed in this paper, only fittings 
deemed "most reliable" were kept. "Most reliable" meaning those which had at 
least 6 or more energy bins for footpoint spectral fitting, and a with a best-fit 

value of 5 or less. Ultimately, 33 "Peak Time" and 89 "Whole Peaks" events 
were used to produce the results that we analyse here. 

To limit pulse pile-up issues (?), care has been taken to discard times with high 
count rates {i.e. just before shutters moving in). Moreover, as spectral fitting was 
usually done above 25 keV, only times with very strong emission (during which 
both attenuators are "in" , or "A3" state) can potentially produce an additional 
component around 35 keV (with only the thin shutter in ("Al" state), detector 
countrates peak around 12 keV, and these can be pile-up to «24 keV. In A3 state, 
the peak of the response is around sslS keV counts. These photons can pile up 
and appear as ~36 keV photons). In 20 of our events were the contribution of 
pile- up photons in certain energy channels larger than 15%. As pile- up typically 
makes two thermal photons appear as a single higher energy photon, imaging 
piled-up photons would place them at the location of the thermal source. In 
19 of these 20 cases, the thermal source was spatially distinct from the HXR 
footpoints, thereby little influencing our results. In the remaining case, the 
thermal source overlapped with the non-thermal HXR footpoints (within our 
7" spatial resolution), and spectral fitting was done above 40 kcV to eliminate 
any contamination by piled-up low-energy photons. 

Table 2 is a list of all parameters obtained for each of our events. Subscripts 
1 and 2 refer to the value of the leading and trailing footpoints, respectively (as 
determined by their heliographic longitude). In a few cases, the subscripts strong 
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Table 2. Measured and derived quantities. 



Symbol Name or decription 



71: 72 


Spectral indices of both footpoints, as obtained by fitting 
a power-law using the OSPEX from the Solarsoft suite of 
routines 


A7 


Spectral index difference A7 = 71—72 between 
footpoints 


7 


= § (71 + 72 ) : average spectral index 


FsO.l, F50,2 


50-keV photon flux in both footpoints [ph/s/cm^/kcV] 


Fso.tot 


=F50,i+F50,2, total flux 


FsO.r 


_ 5(M . 5Q_jjgy fjyx ratio between footpoints 

^50,2 


dt 


Duration or accumulation time [s] 


GOES 


GOES X-ray class, or flux [W m~^] in the low l-SA channel 



Lati, Lat2 Hcliographic longitude [degrees] of both footpoints. 

The footpoint with the largest longitude, (or "leading") is 
labelled "1", the other one "2" 
Loni , Lon2 Hcliographic latitude of both footpoints [degrees] 
s Spherical separation between footpoints [Mm] 

a Angle between footpoints and solar equator 



and weak were also used. They refer to the value of the strongest and weakest 
footpoints, respectively (as determined by their 50-keV flux). 

Presenting all possible combination of scatter plots is prohibitive (they can 
be all found at a website^). Only the most relevant have been presented, but all 
have of course been examined, and an exhaustive table of computed correlation 
coefficients can be found in Section 4.3. 



4. Observations &: Discussion 



4.1. Spatial information 



Flare distribution on the sun: Figure 2 shows the spatial distribution of 
our 53 flares on the Sun. As already known, flares occur predominently at ±15 
degrees of latitude, and there is no marked longitudinal dependence. The slight 
lack of events at high longitudes is very probably due to observational bias: with 
our imaging method (CLEAN with detectors 3 and above), flares with footpoint 
separation smaller than «10" appear to be single-footpoint flare, and are not 
selected. Projection effects near the solar limb reduces the apparent footpoint 
separation, causing some of these flares to be discarded. 
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Flare distribution on the Sun 
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Figure 2. Distribution of our sample of 53 flares on the Sun. 



40 



30 - 




Figure 3. Histogram of spectral indices 7. The gray bars are for "Whole Peaks" events, the 
black ones for "Peak Times" events. The bin size (0.2) was taken to be larger than the average 
error of 0.16. 



FP.tex; 21/11/2011; 1:27; p. 7 



Saint-Hilaire et al. 



Table 3. Consistency/inconsistency of A7 with zero: The numbers on the left of the 
"I" are the number of cases where A7 — n ■ cr^^ >0, and numbers on the right arc 
events where A7 + n ■ a^~f <0. The total number of each type of event is given in 
parenthesis at the top of each column. 



Event type 


"Peak Times" 


"Whole Peaks" 


main "Whole Peak" 




(33 events) 


(87 events) 


of each flare 








(37 events) 


n=0 


17|16 


54|33 


21|16 


n=l, "1-sigma results" 


10|15 


34|19 


12|11 


n=2, "2-sigma results" 


5|4 


16|8 


6|4 


n=3, "3-sigma results" 


2|1 


5|2 


1|0 



4.2. Spectral Information 

Figure 3 shows that spectral indices are generally harder for "Peak Times" 
events, which is of course no surprise, as it is a natural consequence of the Soft- 
Hard-Soft behaviour observed in a majority of flares, where the flattest spectral 
index corresponds to the time of most intense HXR emission (see e.g. ? (?)). 
Flares are very seldom harder than 7 w2.4 in photon spectral index (see e.g. ? 
(?) and references therein). The distribution at high 7 in Figure 3 is not to be 
trusted, as it is distorted by observational bias: only flares with sufficient HXR 
emission above 50 keV , i.e. flares with hard spectra, were used in our study. 

Another observational constraint is the instrument's dynamic range DR: The 
weakest footpoint is visible if F^oak > ^'dr^ • Using a conservative dynamic 
range of «5 for RHESSI, it means that if a footpoint is weaker than the other 
one by a factor 5 or more, it will not be imaged. 

Figure 4 displays scatter plots of the average spectral indices (7) or spectral 
index differences (A7) versus the total 50-keV flux (F5o,tot) or the 50-keV flux 
ratio (Fso.r) of our events, with error bars. And indeed, no event shows a flux 
ratio greater than 5 (or smaller than 0.2). Very few footpoint pairs have spectral 
index difference greater than 0.6, and none above 0.8. This fact could not be 
attributed to observational effects. During our data reduction, a few events with 
A7 larger than 1 were found, but they were discarded because of poor statistics 
(large fltting parameter) and/or the appearance of a third source. 

Roughly 25% of "Peak time" events (8 out of 33) have spectral index dif- 
ferences consistent with zero (i.e., A7 within l-cr of zero; Table 3). This ratio 
increases to sa40% for "Whole Peak" events (34/87). 

No "Peak time" event (0/33, Table 4) has a flux ratio consistent with unity, 
and only salO% of the "Whole peak" events do (8/87; this fraction is almost 
reduced to zero (1/37) when considering the strongest "Whole peak" events of 
each flare {i.e. those encompassing the "Peak time")). Table 4 seem to suggest 
that leading footpoints might have more flux during "Peak time" events, but the 
result is not statistically signiflcant, and will not be further discussed. 



http: / /sprg. ssl.berkeley.edu/~shilaire/FootPointProject/htmlsummaries/browser. html 
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Figure 4. Top: 7 vs. F50, tot, with 1-cr error bars. Bottom: A7 and IA7I vs. Fso.r, with l-cr 
error bars (see Table 2 for an explanation of all quantities). The data points in gray are the 
ones for which | A7I < cr/^^. The thick, gray line is a linear regression to the A7 vs. Fso^r "Peak 
times" data. 



Figure 5 and Table 3 show that there is no statistically significant preference 
for the leading footpoint to be either harder or softer than the trailing one, 
during either "Peak Times" or "Whole Peaks" . 

The deficit of F.so.r ~1 events during "Peak Times" is particularly clear in 
Figures 4 (bottom left) and 6 (only 1 out of 33 events is within 10% of unity 
flux ratio, and only 5 out of 33 within 20%). This greater footpoint asymmetry 
during times of peak HXR fluxes than during whole HXR peaks could indicate 
that individual particle acceleration episode occur preferentially in one direction 
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-1.0 -0.5 0,0 0.5 1.0 0.0 0,2 0.^ 0,6 0,8 1,0 
A-/ = 7,-72 IA7I = 17,-72 



Figure 5. Histogram of spectral index differenees A7. Crosshatciied: "Peak Time" events. 
Solid black: "Whole Peak" events. Solid gray: "Whole Peak" events (only one per flare, the 
one overlapping the peak HXR flux time). The bin size (0.3) was taken larger than the average 
error (0.23). 




50,1/^50,2 ^50,slrong/^50,weok 



Figure 6. Histogram of flux ratios. Crosshatshed: "Peak Time" events. Solid black: "Whole 
Peak" events. Solid gray: "Whole Peak" events (only one per flare, the one overlapping the 
peak HXR flux time). The bin size (0.5) is much larger than the typical error. 
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Table 4. Consistency /inconsistency of with unity: Events with Fj — n ■ crp^ >1 are 
left of the "I", and events with Fr + n ■ crp^ <1 are at the right. 



Event type 


"Peak Times" 


"Whole Peaks" 


main "Whole Peak" 




(33 events) 


(87 events) 


of each flare 








(37 events) 


n=0 


20|13 


41|46 


21|16 


n=l, "1-sigma results" 


20|13 


40|39 


21|15 


n=2, "2-sigma results" 


19|12 


32|34 


16|13 


n=3, "3-sigma results" 


19|11 


24|32 


12|12 



Peak Times" "Whole Peaks 




1 2 3 4 5 1 2 3 4 5 

^50, strong/ ^50. weak ^50, strong/ ^50, weak 

Figure 7. F5o,strong/-F50,weak VS. 7strong - 7wcak , where strong (weak) denotes the footpoint 
with the strongest (weakest) 50-keV flux, respectively. The dashed line marks zero spectral 
index difference. There are 23 out of 33 (79%) "Peak flux" events and 53 out of 87 (61%) 
"Whole peak" events which lie below the dashed line (see also Table 5.). 

of the loop at any given time, but that, on the average, particles tend to be 
accelerated in both directions equally. 

Table 5. Column density model; (agreement/disagreement). 



Event type 


"Peak Times" 


"Whole Peaks" 


All together 


All events 


30.3% (10/33) 


39.1% (34/87) 


36.7% (44/120) 


1-sigma results 


24% (6/25) 


42.9% (21/49) 


36.5% (27/74) 


2-sigma results 


28.6% (2/7) 


40.0% (8/20) 


37.0% (10/27) 


3-sigma results 


33% (1/3) 


60.0% (3/5) 


50.0% (4/8) 
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We checked in a simple way the agreement of our data set with the theory 
presented in Appendix A: that spectral index differences A7 between flare foot- 
points might be due to differences in column depths in asymmetric loops. This 
effect, assuming equal distributions of electrons are accelerated in both directions 
of the loop, results in the footpoint having the softer spectrum also having the 
most flux. As shown in Figure 7 and Table 5, the reverse happens most of the 
time, i.e. better supporting a "magnetic mirroring" type of effect (see e.g. recent 
work by ?, ?, and references therein). Furthermore, the best candidate in support 
of the column density difference model is the 23 July 2002 flare (which agrees at 
the "3-cr level") actually hits a snag when one considers the amount of 50-kcV 
flux (predicted by theory) out of the footpoints (i.e. emitted somewhere along 
the legs of the loop, before reaching the imaged footpoints): as explained in 
Appendix A, this 50-keV emission should have been observable. Moreover, the 
presence of large flux ratios F^ (2 or above, ^ or below) also lead us to believe 
that this theory cannot be a dominant factor, at least for reasonable values of 
leg column densities (see Appendix A). 

There are several altenatives to explain footpoint asymmetries: (a) Asym- 
metrical acceleration: The strong footpoint asymmetries, particularly during 
HXR peak times, suggest that it is the acceleration mechanism itself which 
could be asymmetrical. If the acceleration process actually took place in the 
chromospheric footpoints (as opposed to high in the corona), one would expect 
asymmetries, as both acceleration processes could in principle be independent 
from one another, (b) Non-uniform target ionization: It is conceivable that both 
chromospheric footpoints have different ionization structure, i.e. considering the 
simplified step model of ? (?) or ? (?), that the column density required before 
reaching the lower-chromospheric regions of neutral atoms is different in both 
legs of the loop, perhaps due to some prior heating of only one of the footpoints. 
Modeling and comparison with observations are required to further this idea, 
(c) The best candidate mechanism to explain footpoint emission asymmetries 
is magnetic mirroring, as discussed in ? (?)'s trap-hprecipitation model (?): the 
magnetic field converges more rapidly in one of the footpoints, and particles are 
mirrored before they reach the dense lower regions (see also ?, ?). Our data 
corroborate better that scenario than the column density asymmetry model. 
The effects of photospheric albedo (?) or return currents (?) might reinforce any 
asymmetry observed in footpoint photon spectra, but only if the accelerated 
electron distributions started out different. 

4.3. Correlation Table: 

Correlation coefficients have been computed for our data, and arc displayed in 
Table 6. 

We found the following: 

• There is some degree of anti-correlation (-0.53, with the 95% confidence 
interval being [-0.74,-0.22]) between A7 and Fr, for "Peak Times" events. 
This is a consequence of there having more events in the upper left and 
lower right quadrants of the lower left panel of Figure 4, and has already 
been discussed in Section 4.2. 
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flux ratio 

^50.2 I 
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ln{F50,l + -F50,2) 



A7 



IA7I 
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Duration [s] 
' f50.l \ 
,^50,2 J 
ln(-F50~i + -F50,2) 
A7 

IA7I 
7 

Longitude 
I Longitude I 
Latitude 
|Latitude| 
FP separation 

FP angle 
In(GOESflux) 



-/1. 00 



-/-0.08 
1.00/1.00 



-/-O.Ol 
0.06/0.12 
1.00/1.00 



-/-0.07 

-0.53/-0.23 

0.24/0.11 
1.00/1.00 



-/-O.ll 

0.45/0.09 

-0.22/-0.13 
-0.42/0.16 
1.00/1.00 



-/-0.21 

0.32/0.22 

-0.23/-0.24 
-0.41/-0.08 
O.5O/0.49 
1.00/1.00 













FP 


FP 


GOBS 




Longitude 


Longitude 


Latitude 


1 Latitude 


separation 


angle 


class 


Duration [s] 


-/0.02 


-/-0.08 


-/-0.03 


-/O.IO 


-/0.52 


-/-0.05 


-/0.22 


V ^50,2 J 


0.15/0.05 


-0.24/-0.12 


-0.01/0.07 


0.43/0.10 


-0.12/-0.07 


-0.07/0.01 


0.16/0.12 




0.26/-0.06 


0.12/0.17 


0.20/0.01 


-0.25/-0.14 


0.04/0.00 


0.22/0.20 


0.88/0.71 


A7 


-0.01/-0.08 


0.33/0.31 


-0.20/-0.19 


-0.12/-0.18 


-0.18/-0.03 


0.05/0.25 


0.11/0.03 


IA7I 


-0.01/-0.02 


-0.25/0.12 


0.12/-0.20 


0.19/0.13 


0.01/-0.18 


-0.04/0.08 


-0.09/-0.07 


7 


0.29/0.31 


-0.24/-0.16 


0.06/-0.03 


-0.06/0.11 


-0.01/-0.42 


-0.26/-0.02 


-0.12/-0.22 


Longitude 


1.00 


-0.04 


0.17 


-0.06 


-0.02 


0.09 


0.12 


1 Longitude] 




1.00 
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Figure 8. "Whole Peak" duration vs. footpoint separation. The error bars are small and were 
not omitted. The gray line is a linear fit (using the bisector method) to the data. It yields a 
power-law slope of O.GitO.l. 



• 7 and IA7I seem also slightly correlated: indicating that IA7I is larger when 
the flare is softer. Upon closer examination, it appears softer spectra to be 
a simple case of softer spectra having larger errors. 

• An unexpected correlation -albeit weak (the 95% confidence interval is 
[0.32, 0.67])- was fomid between footpoint separation and event duration, 
for "Whole Peaks" events. 

From Figure 8, it seems that we have 

(HXR burst duration) « (footpoint separation)^, (1) 
or, assuming semi-circular loops: 

(HXR burst duration) « (loop length)^ (2) 

It seems that the longer the loops, the longer the HXR peak will last. The 
interpretation is not yet clear. It could be a simple case of larger loops 
needing more time to evolve than the short ones during the flaring process, 
or another case of the "big flare syndrome" : everything is bigger in larger 
flares. 

• Excellent correlation between "Peak Times" total flux and GOES class 
(95% confidence interval for correlation coefficient: 0.77-0.94), less good 
for "Whole Peaks" total flux and GOES class (95% confidence interval for 
correlation coefficient: 0.57-0.81). 

Figure 9 shows the clear correlation between event's maximum GOES 1-8A 
flux and the total HXR flux F50, for "Peak Time" events. Fitting a power-law 
yields: 
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Peak Times 



1 00.0 



^ E 

+ CJ 



10.0 



1 .0 



0.1 I . . . 

10"^ 10"^ 10"'- lO"-' 

GOES max flux [W/m^] 

Figure 9. GOES maximum flux in the l-SA band vs. peak 50-keV flux of both footpoints, 
with 1-tT errors. The solid line is a power-law linear fit to the data. The dotted line is another 
power-law fit to the data, using the bisector method. 



F50 = A • F" o (3) 

GOES,l-8A 

where A = (4.7 ± 0.3) x 10^ and a = 0.8 ± 0.1 when F50 is in photons s~^ cm~^ 
keV~^ and F ° in W m~^. The bisector method (?) is more relevant 

GOES,l-8A ^ ' 

when the variables are independent, in which case we have A = io*-3±0-2 a^d 
a = 0.97 ± 0.05. (Of course, "Whole HXR Peaks" events typically lie below this 
solid line, with a wide scatter.) 

This good correlation can be interpreted as a direct consequence of the Ne- 
upert effect (?; ?): the larger the amount of non-thermal energy (approximated 
by F50), the larger the amount of thermal energy (approximated by the GOES 
SXR) (since F50 is a power-law normalization factor, it describes equally well 
the amount of non-thermal electrons of lower energies, which contain most of 
the non-thermal power). Similar correlations have been reported before (see e.g. 
?, ? and references therein.) 



5. Summary and Conclusion 

The following is a compilation of our results, and can be used as a list of 
contraints for any flare and particle acceleration theories: 

(1) The total footpoint 50-keV flux correlates remarkably well with the GOES 
maximum l-SA flux. The relationship is fairly linear. 

(2) There is no statistically significant difference in our sample between "lead- 
ing" and "trailing" footpoints, as regards asymmetries. 
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(3) Flares are mostly located at ±15 degrees of solar latitude, and flare 
parameters have no marked longitudinal dependence. 

(4) Spherical separation between footpoints seem not to correlate with any 
of the other parameters examined, with the surprising exception of HXR burst 
duration, where a weak correlation was found. This seems to indicate that longer 
loops produce longer HXR peaks, probably because the magnetic disturbance 
and particle acceleration last longer in long loops than in short ones. 

(5) Flare footpoint spectral indices 7 are seldom below «2.4 (1 case out of 
172). "Peak times" are generally harder than "Whole Peak" intervals, a natural 
consequence of the commonly observed soft-hard-soft behaviour of flares. 

(6) w25% ("Peak times") to sa40% ("Whole Peaks" ) of double footpoint flares 
have spectral index differences A7 consistent with zero. A7 can reach 0.6, and 
only rarely goes beyond. The amplitude of A7 is uncorrelated with flare GOES 
class. 

(7) 50-keV footpoint flux ratios are never quite unity, are typically between 1 
and 2, and only seldom go beyond 3. This result could be due to observational 
bias. 

(8) The asymmetric loop model, where a column density difference is responsi- 
ble for the difference in spectral index and flux between HXR footpoints, cannot 
explain a majority of our observations. It is therefore not a dominant factor. 

(9) The greatest asymmetry being around "Peak Times" further suggests that 
magnetic reconflguration is greatest at those times. 

Appendix 

A. Column Density Effects in Asymmetric Loops 

Figures 10 and 11 show the numerically-computed effect of loop asymmetry on 
the 50-kcV component of the thick-target brcmsstrahlung spectra produced by 
two identical accelerated electron distributions. 

Notice that emission at the footpoint of the leg of the loop with the least 
column density is softer than the other one, but has more flux. 

Increasing S increases all three parameters (A7, Fso^r, and Coronal/FP flux), 
whereas increasing / only increases the first two, and makes the last one decrease. 
Hence, in order to have large A7 and for the spatially extended coronal flux to be 
lost within the dynamic range of the instrument, a high / (i.e loop asymmetry) 
is required. 

Figure 12 show the same data as Figures. 10 & 11, but with different axes. 
The bottom plot of Figure 12 clearly shows that, in order for a hard flare 
(such as we have in this study) to have a flux ratio beyond 2 (or smaller than 
1/2), and have the coronal part of the thick-target 50-keV flux go unobserved, 
unreasonable values for / (beyond 10!) or N2 would have to be considered. In 
the trap-|-precipitation model, the footpoint with the highest (in altitude) mirror 
point will stop only the lowest-energy electrons. All the hight-energy electrons 
will mostly stop (and emit brcmsstrahlung) in the other footpoint. In this case, 
the footpoint with the hardest emission will also have the most flux. 
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Figure 10. If identical accelerated electron distributions travel through two different col- 
umn densities A^l and N2 (such as would be the case in an asymmetric loop), then the 
bremsstrahlung photon spectra they will emit at each footpoint will be slightly different. Top: 
Spectral index difference A7 = 71—72 between footpoint spectra, assuming N2=2xNl. 



Middle: 50-kcV flux ratio F50 r 



Fso.l 
F50.; 



between footpoints. Bottom: Ratio of the total 50-keV 



flux present in the corona to that present in both footpoints. The footpoint spectral indices 
and 50-keV fluxes were determined by fltting a power-law in the 25-75 keV band, typical of the 
observations presented in this paper. The hatched part indicates areas where 50-keV coronal 
emission would be observed, assuming a conservative dynamic range of 5 for the instrument. 
5 is the accelerated electron power-law spectral index. 



B. RHESSI Imaging Spectroscopy Errors 

Errors in RHESSI imaging spectroscopy are extrememly difficult to estimate, as 
each individual pixel or feature in an image is heavily correlated to other parts 
of the image, via the point-spread function. 

The current heuristic method implemented in OSPEX, the standard RHESSI 
imaging spectroscopy software package, is to define the error on the flux of a 
feature as the maximum of the flux outside of the (visually-selected) sources, 
divided by a somewhat arbitrary value of n =3. When the data is very noisy 
{e.g. images made at high energies, were the counts are low), this method will 
assuredly under-estimate the errors. This later point is not an issue for our study, 
though, as we have discarted noisy images with our choice of energy bands. 

We have tried another method, similar to Mitani (2005), where a gaussian 
is fitted to a histogram of the pixel values of the whole image. The error on a 
source flux is then taken to be the l-tr extend of this distribution, multiplied by 
the number of pixels in that source. We have labelled that method "Gaussian 
Background" . 
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To estimate the accuracy of each methods' error estimation, we have plotted 
(sec Figure 13) the histogram of the normahzcd residuals of all our spectral 
fittings. Assuming our power-law model is correct, then the leftmost plot of 
Figure 13 leads us to conclude that errors in the standard OSPEX imaging 
spectroscopy package arc typically over-estimated by a factor w2. Wc then chose 
n =6 as our heuristic number, and obtained the plot in the middle of Figure 13, 
where a fitted gaussian has a very close to unity. This is the scheme that we 
finally settled upon. 

The fitting parameters hardly change at all whether we choose to use OSPEX 
with n = 3 or 6, or with the Gaussian background method, only the l-cr error 
on those parameters are infiuenced by the choice of the method. Moreover, the 
results and conclusions obtained using either the OSPEX with n =6 or the 
Gaussian backgound method do not change significantly. 
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Normolized residuols Normalized residuols Normolized residuals 

Figure 13. Histogram of normalized residuals. Left: OSPEX with n =3, middle: OSPEX with 
n =6, and right: Gaussian background method. is the total number of data points fitted, fj, 
and cr the mean and standard deviation of the gaussians fitted to the histograms. 
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